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ABSTRACT 

We consider the gravitational stability of a current carrying filamentary cloud in the pres- 
ence of both axial and azimuthal magnetic fields using a simple analytic model. The azimuthal 
magnetic field is shown to give rise to a new contribution, dictated by Ampere's law, in the 
corresponding virial equation for magnetohydrodynamic equilibrium. From this we obtain a 
computationally inexpensive guidance on the gravitational stability of current carrying filamen- 
tary clouds. The approach not only provides a fresh insight into the essential physical mechanisms 
involved, but also demonstrates clearly that, for sufficiently large and yet astronomically realistic 
currents, the azimuthal magnetic field can cause filamentary clouds to undergo instability. 

Subject headings: ISM: clouds - ISM: magnetic fields - ISM: structure - stars: formation 



1. INTRODUCTION 

One of the fundamental conditions for the onset 
of the dynamical instability causing the disperse, 
cold gas to collapse was given by the Jeans con- 
dition (Jeans 1902). This gave us a requirement 
for the instability to occur as an expression of the 
velocity of a propagating wave through the gas. 
When the velocity of the wave became imaginary 
we would no longer see oscillations but an expo- 
nential collapse at the boundary instead. This 
could also be expressed as a minimum, critical 
wavelength beyond which we would see an insta- 
bility occurring. When considering these molec- 
ular clouds they also arise in cylindrical filamen- 
tary mass distributions, such as the "Nessie" Neb- 
ula (Jackson et al. 2010). Further, infrared dark 
clouds (IRDCs) with a filamentary structure have 
been identified by new telescopes such as Herschel 
(Arzoumanian et al. 2011; Hill et al. 2011). 

To account for this matter distribution the 
Jeans stability condition was modified. This mod- 
ification showed that cylindrical clouds, subject 
to longitudinal perturbations, can (under specific 



conditions) break apart into smaller pieces (Dibai 
1957), the results of which could be used to explain 
the creation of star chains. Both of these methods 
to understand this instability neglect the magnetic 
field and its influence on the stability of the sys- 
tem. It was not until 1953 when Chandrasekhar 
and Fermi approached this problem for the first 
time (Chandrasekhar & Fermi 1953). They re- 
alized that for stationary configurations it is suf- 
ficient to consider a less stringent stability con- 
dition called the virial equation. The method of 
analyzing the stability through the virial equation 
has the advantage of being much simpler than di- 
rectly tackling often cumbersome partial differen- 
tial equations, this is often the case when consid- 
ering a method similar to the derivation of the 
Jeans condition (Jeans 1928). By considering a 
plasma, they first derived the virial equation in 
the most general case including a magnetic en- 
ergy term. This was done by first considering the 
Navier-Stokes equation in three dimensions with 
the addition of magnetic force density terms. By 
integrating the equation over the volume they ob- 
tain an expression involving the energy terms from 



the different components, i.e, the magnetic energy 
and the molecular energy. This can be used to rep- 
resent a condition required for the system to be in 
equilibrium. Using a similarly fashioned deriva- 
tion as done in the most general case, they then 
did the same for cylindrical distributions of matter 
with a prevailing axial magnetic field. From this 
a stability condition was reached. It was shown 
that the average magnetic intensity is required to 
be less than a critical value, which is dependent 
on the average mass density of the cylinder, in or- 
der to maintain stability. However, Dibai (1958) 
showed that for longitudinal waves propagating 
through a cylindrical medium of plasma, an ax- 
ial magnetic field does not affect the stability of 
the system. 

This previous research focuses on analytical 
methods by first assuming that the system starts 
close to gravitational stability and then showing it 
can fragment. However, the interstellar medium 
(ISM) is highly dynamical and yet filaments can 
be shown to form in these environments through 
the use of numerical simulations (Ntormousi et 
al. 2011). The gravitational collapse in turbu- 
lent molecular clouds is influenced by the magnetic 
field. It can be shown that magnetic fields cannot 
prevent a local collapse unless they act as a mag- 
netostatic support (Heitsch ct al. 2001). Despite 
the turbulence of the ISM it is possible to form 
pillar like molecular clouds and through numerical 
simulations it is shown that it can be due to im- 
pact ionization from nearby massive stars. Grav- 
itational instability can be induced at the tip of 
these pillars leading to the formation of low mass 
stars (Gritschneder et al. 2009). Molecular clouds 
in general can form in many different shapes and 
sizes but the background magnetic field can be 
shown to strongly influence its development. Fur- 
ther numerical simulations show that the direction 
and strength of this magnetic field is vital in deter- 
mining the form of the molecular cloud (Heitsch 
et al. 2008). The properties of clumps of interstel- 
lar material formed in the presence of weak mag- 
netic fields with a colliding warm neutral medium 
is investigated through the use of numerical mag- 
nctohydrodynamic (MHD) simulations (Banerjce 
et al. 2009). It was shown that the magnetic field 
has an alignment with the velocity but there exist 
intense fluctuations both in magnitude and direc- 
tion within the clumps meaning that the magnetic 



field has become distorted by the turbulence of the 
gas. Despite this a dynamical collapse of the cloud 
is still possible. The evolution of giant molecular 
clouds when subject to a collision with a warm 
neutral medium streams in the presence of mag- 
netic fields and ambipolar diffusion has been inves- 
tigated through the use of numerical MHD simu- 
lations (Vazquez et al. 2011). It was shown that 
only for supercritical or marginal critical streams 
could reasonable rates of star formation occur. It 
was also shown that the mass to magnetic flux ra- 
tio is a fluctuating function of position within the 
cloud which they use to imply that a large por- 
tion of the cloud's mass may remain magnetically 
supported against collapse whereas star formation 
occurs in the non-supported regions. 

There have been observations to suggest that 
some molecular clouds may exhibit helical mag- 
netic fields (Bally 1987), consisting both axial 
(poloidal) and azimuthal (toroidal) components. 
It has also been found that a current carrying jet 
creating an azimuthal magnetic field is a useful 
model to explain the temporary capture of car- 
bon monoxide (CO) in the spiral galaxy NGC4258 
(Krause et al 2007). The presence of the az- 
imuthal magnetic field will influence the stabil- 
ity of the system, which has been addressed by a 
number of recent papers (Fiege & Pudritz 2000a, 
2000b). The effect of the azimuthal magnetic field 
on rotating, isothermal clouds has been investi- 
gated in (Tomisaka 1991). The main results of 
these works are obtained through MHD numerical 
simulations. The gravitational collapse and frag- 
mentation of magnetized filamentary clouds has 
been investigated and it can be shown that an az- 
imuthal magnetic field can strengthen the cylin- 
drical cloud against fragmentation during collapse 
(Tilley & Pudritz 2003). Numerical MHD simu- 
lations including turbulence have shown that the 
magnetic field plays a significant role in the col- 
lapse of molecular clouds and the formation of 
star-forming cores (Tilley & Pudritz 2007). The 
influence of the magnetic field on the star forma- 
tion rate in molecular clouds is profound. Even in 
the presence of weak magnetic fields the formation 
of dense cores in the star formation process can be 
slowed significantly resulting in a longer evolution 
time (Nakamura & Li 2011). 

The purpose of the present paper is to provide 
an analytical approach to the gravitational sta- 



2 



bility of a cylindrical filament of plasma which 
is subject to both axial and azimuthal magnetic 
fields. We focus on the essential physical mecha- 
nism and aim to provide a computationally inex- 
pensive guidance on the corresponding necessary 
condition of gravitational stability. Starting from 
the most general fluid flow equation we derive a 
new form of virial equation. Specifically, it yields 
information on how gravitational stability is af- 
fected by the presence of azimuthal magnetic field. 

2. THE CONFIGURATION OF A CYLIN- 
DRICAL PLASMA WITH MAGNETIC 
FIELDS 

We shall idealize a dark filamentary plasma as 
an infinitely long cylinder with radius R of fluid 
having zero resistivity, along whose axis we allow 
for a net movement of charge. This sets up an 
azimuthal magnetic field, just like that of a cur- 
rent carrying wire, in addition to a possible pre- 
vailing magnetic field in the axial direction. It is 
natural to adopt the cylindrical polar coordinates 
(r, <j) 1 z) where z denotes the axis of the filament. 
The cylindrical symmetry requires the all physical 
quantities involved including the density p, pres- 
sure p, velocity v and magnetic field B to be func- 
tions of radial distance from the axis r and time 
t. In particular both azimuthal and axial com- 
ponents of the magnetic field, B$ and B z respec- 
tively, depend only on r and t. 

The magnetic fields across the boundary of the 
cylinder are assumed to match smoothly (contin- 
uous and diffcrentiablc) so that there is no surface 
current sheet. Exterior to the cylinder the axial 
magnetic field B z is assumed to be zero for sim- 
plicity while the azimuthal magnetic field B^ may 
remain non-vanishing. 

3. VIRIAL EQUATION WITH AXIAL 
AND AZIMUTHAL MAGNETIC FIELDS 

We begin by considering the Navier-Stokes 
equation for an inviscid fluid in the MHD ap- 
proximation as follows 
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where v is the velocity, t is the time, p is the pres- 
sure, V is the gravitational potential, and is the 
p density . Here 4? denotes the Lagrangian time 



derivative of v. The cylindrical symmetry to be 
used throughout yields 
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and 
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2p, or 2p 

in cylindrical polar coordinates. The second term 
on the right-hand side of Equation ([lj is a gravi- 
tational term, satisfying 
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where g is the gravitational acceleration vector 
with radial only components at radius r , dA is the 
surface area unit, and M onc is the mass enclosed 
by an imaginary cylindrical surface of radius r and 
length L. This yields 
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where A onc = M cnc /L is the mass per unit length 
enclosed within radius r. For r > R we shall de- 
note the constant mass per unit length by A = 
Acne- Thus we see that 
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The third term on the right-hand side of Equa- 
tion ([1]) is the magnetic tension term. To evaluate 
this term, we consider the following general expan- 
sion in cylindrical coordinates: 

(P.V)Q = 
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for two arbitrary vector fields P and Q. On re- 
stricting both P and Q to be B = Bfi(r,t)(j) + 
B z (r,t)z, the above expansion reduces simply to 
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See, e.g., Arfken (1985). 

Therefore, by substituting relations ©, ©, © 
and (O into Equation (TTJ), we see that 

dv r dv z d |B| 2 

2GA C nc » ^6 ~ /0 n 

pr — r. (8) 
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Integrating Equation ((SJ over the cross-section 
area of the cylinder, we obtain two equations from 
the resulting z and f components respectively as 
follows: 
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To proceed we shall denote by 
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where dX is a mass per unit length element. Eval- 
uating the gravitational integral first, we get 
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We evaluate the first integral on the right-hand 
side of Equation (fTOj) by considering the angular 
integral, 
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We then perform the radial integration by using 
the method of integration by, parts, 

-2. r R r 2 (^){ P +^)dr = 
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So the first integral on the right-hand side of 
Equation (fTU|) becomes, after integration, 
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where p b is the pressure at the boundary, which is 
zero providing we have no surface currents. This 
is justified with the smooth matching of the mag- 
netic field at the boundary. 

Together with the simplifying assumption B z = 
at the boundary, we can rewrite Equation (|15[) 
in the following form: 
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where U is the thermal energy of motion per unit 
length, 7 is the ratio of specific heats, and M. is 
the total magnetic energy per unit length. The 
subscripts on M. indicate the contributions from 
the magnetic field components, i.e., azimuthal and 
axial. 

To carry out the last integral on the right-hand 
side of Equation (flOl) , 
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we substitute Equations ([12]) , (TT6l) and (fT7f into 
Equation (|10j) and consider equilibrium with the 
acceleration term equal to zero. This leads, after 
some algebra, to the equation 

B 2 
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where is evaluated on the surface of the cylin- 
der at r = R, and 

Q = -GX 2 (19) 

denotes the gravitational energy per unit length of 
the cylindrical cloud. Interestingly, the azimuthal 
magnetic energy M.$ due to the magnetic tension 
in Equation (|17l) is canceled by a similar term with 
a minus sign due to the magnetic pressure in Equa- 
tion erg) . 

By using Ampere's law 

s H^) 2 < 2 °> 

in terms of the total current / carried by the cylin- 
der, we see that Equation (|18p becomes 

Ct + 2( 1 -l)U + 2M z -^I 2 =0. (21) 

An 

Equation (|21j) is our new virial equation, repre- 
senting hydrostatic equilibrium in the presence of 
both axial and azimuthal magnetic fields. 

It is worth noting that if we set the current to 
zero and remove the azimuthal magnetic field then 
we recover the virial equation 

+ 2(7 - l)U + 2M z = (22) 

due originally to Chandrasekhar & Fermi (1953) 
for the case with only axial magnetic field. 

Remarkably, the last term — ^ I 2 of Equation 
(f2~Tj) featured in this new virial equation, arises 
from the boundary term when integrating the 
magnetic pressure. Indeed, the existence of this 
new contribution is dictated by Ampere's law for 
a current carrying filamentary cloud, and is absent 
from a finite three-dimensional cloud. Further- 
more, its minus sign has an important implication 
that while the magnetic field tends to stabilize a fi- 
nite three-dimensional cloud, see, e.g., Shu (1992), 
the azimuthal magnetic field can destabilize a fil- 
amentary cloud. This is to be expected through 
the pinch effect of MHD, which we shall explore 
in the next section. 

4. GRAVITATIONAL STABILITY WITH 
AZIMUTHAL MAGNETIC FIELD 
ONLY 

Here we shall set the axial magnetic field to 
be zero and focus on the effect of the azimuthal 



magnetic field. Thus, 

M z =0. (23) 

The aim is to provide a simple estimate of the 
gravitational stability of the system based on the 
strength of the azimuthal magnetic field and mass 
density of the cloud. Therefore, we consider an 
idealized case where the mass density p and cur- 
rent density J are homogeneous throughout the 
medium. 

The total energy per unit length of the system 
takes the form 

E = n + U + M4, (24) 

where Aitf, is the azimuthal magnetic energy per 
unit length. In our present estimate, we shall eval- 
uate M.cf, using the azimuthal magnetic field inte- 
rior to the cylinder: 

M+= / —rd^dr. (25) 

J0=O Jr=0 

by assuming that the azimuthal magnetic field is 
effectively canceled by that of the return current 
exterior to the cylinder (Falgarone & Passot 2003). 

The integral in Equation (|25|) can be readily 
evaluated using the constant current density J and 
the corresponding total current I = irR 2 J, yield- 
ing 

M,-^RJ - r^I ■ (26) 

Therefore, by applying Equations (|2"3"]l and (|2lfl l , 
the virial equation (|2i"j) now becomes 

fi-4M* + 2(7 -1)17 = 0. (27) 

Eliminating U from Equation (I27|) by using the 
energy expression ([24| we see that 

2( 7 -l)£: = 2(7 + l)7W + (27-3)O. (28) 

As a necessary condition for the dynamical sta- 
bility of the cylindrical cloud, its total energy per 
unit length must be less than zero (E < 0). From 
Equation (|28[) . this means the following inequality: 

-(2 7 -3)fi > 2{ 1 + l)M 4> . (29) 

Therefore, from Equation (|2"91 we can reach a 
necessary condition for stable equilibrium, using, 
e.g., 7 = | for monatomic plasma, as follows: 

-Jr!>AV (30) 
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This places an upper limit on the azimuthal mag- 
netic energy to be a quarter of the absolute gravi- 
tational energy for the system to be gravitationally 
stable. 

By using Equations (fl9l) . (|26|) and the expres- 
sion for mass per unit length A = irR 2 p the in- 
equality (|30|) becomes 

\ G (^ P r > (3D 

leading to the conclusion that for gravitational sta- 
bility the current density must necessarily satisfy 
the following inequality: 



5. DISCUSSIONS 

The new form of virial equation (|21l) gives us a 
gravitationally bound condition for the filament in 
hydrostatic equilibrium, this alone however may 
not mean that the filament is stable. The re- 
sult can be reduced to the one derived by Chan- 
drasekhar & Fermi (1953) by setting the current 
(and hence the azimuthal magnetic field) to zero. 
This is an important result because it demon- 
strates that the azimuthal magnetic field signifi- 
cantly modifies the previously established neces- 
sary condition for gravitational stability. How- 
ever, in making this comparison we need to be cau- 
tious in comparing our result with their findings. 
The azimuthal magnetic field can extend to infin- 
ity whereas the axial magnetic field is restricted 
to existing solely inside the plasma and falling to 
zero at the boundary but it need not be subject to 
such a restriction. It is possible for the axial mag- 
netic field to exist outside the plasma and to be 
non-zero at the boundary. The assumption which 
is made that the pressure of the plasma is zero 
at the boundary need not be the case if surface 
currents exist. This would create a Lorentz force 
at the surface inducing pressure at the boundary. 
These additional considerations will no doubt af- 
fect the gravitational stability of the system and 
should be accounted for. On the other hand even 
without these considerations we can now say with- 
out doubt that an azimuthal magnetic field does 
indeed influence the gravitational stability of the 
system, proportionally to the current squared. 



Removing the axial magnetic field, leaving the 
azimuthal magnetic field, we derive the gravita- 
tional stability condition given by Equation (|32p . 
This inequality (f3"2"j) sets up an estimated upper 
limit of the current density beyond which gravita- 
tional instability occurs, for given density of the 
cylindrical filament. For example, using a rela- 
tively low density of 10 4 cm~ 3 at 10 K (Henning 
et al 2010) and the radius of 0.15 lt-yr (Arzou- 
manian et al. 2011), whose values are typical of 
such filaments, we see that under these conditions, 
inequality (f3"2"j) yields: 

J < 4.4 x 10~ 19 Am- 2 . (33) 

In this case, it is clear that the current density on 
the right-hand side of the above inequality physi- 
cally represents a few electrons per second passing 
through a meter squared of area. Using the given 
radius, this value is equivalent to the minimum 
current of 2.8 x 10 12 A to violate gravitational sta- 
bility. 

The above current generates an azimuthal mag- 
netic field of 3.9/zG, at the boundary of the cylin- 
der, which is around the expected strength of mag- 
netic field. For example, the local ISM has a mag- 
netic field strength of 4/xG (Opher et al. 2009). 
The Radio Arc filaments have been estimated to 
have a magnetic field strength of 10/iG (Ferriere 
2009). Let us examine the Serpens South Clus- 
ter, which is 848 lt-yr from Earth. It is in a very 
early part of the star formation process. It is 
emitting near-infrared light detected by Herschel 
(Bontemps et al. 2010) and NASA's Spitzer Space 
Telescope (Gutermuth et al. 2008). 

This IRDC is predicted to have a magnetic 
field of a few hundred micro Gauss in two dif- 
ferent zones along the main filament (Sugitani et 
al 2011). According to our analysis, it would be 
gravitationally unstable. Consequently, small per- 
turbations around its MHD equilibrium may be 
subject to dynamical instability, thereby growing 
in an unbound manner, while the stable parts of 
the filament remain as fragments. This fragmenta- 
tion could be followed by the eventual collapse of 
the resulting fragments, trigging star formation, in 
fashion of Jeans' instability. Such processes have 
been demonstrated through MHD numerical sim- 
ulations (Fiege & Pudritz 2000a,2000b) and semi- 
analytical solutions for self-similar reduced models 
of MHD (Hennebelle 2003), where certain pertur- 
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bative oscillations of wavelength greater than some 
scale related to the azimuthal magnetic field have 
been observed. This scenario could be relevant 
for the observation of the above IRDC, where the 
infra-red detectors show, at its core, a bright clus- 
ter with more stars forming in a dispersed manner 
along the body of the cylinder. There are 35 pro- 
tostars in total beginning to form (Bontemps et 
al. 2010). 

The idealistic magnetic fields which we have 
considered can be expanded to consider a more 
complex magnetic field geometry, perhaps as a re- 
sult of turbulence (Federrath et al. 2011), it can be 
shown that weak magnetic fields can be exponen- 
tially amplified through small dynamo effects dur- 
ing the process of star formation (Sur et al 2010). 
Such structures producing tangled magnetic fields 
can affect stability, e.g. by making the molecular 
cloud more difficult to compress. 

We have adopted throughout the cylinder to 
be infinitely long to provide a first approximation 
to filamentary clouds with large aspect (length to 
width) ratio. However, it is also possible to extend 
this model to a closed filamentary structure. This 
is where we have a finite radius but a feedback 
mechanism for the current within the filament. In 
this case, there would exist the current inside the 
molecular cloud producing the azimuthal magnetic 
field but in addition there would exist a return cur- 
rent outside of the cloud. The form of this return 
current would depend on whether the molecular 
cloud was isolated or existed in a much larger net- 
work, both of these possibilities have been briefly 
explored (Falgarone & Passot 2003). If a finite 
cylinder is considered higher modes will be ex- 
cited, requiring a more in-depth analysis. 
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